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ABSTRACT 

We argue that the outbursts of the FU Orionis stars occur on timescales which are 
much longer than expected from the standard disc instability model with a c > 10 -3 . 
The outburst, recurrence, and rise times are consistent with the idea that the accretion 
disc in these objects is truncated at a radius Ri ~ 40 Kq. In agreement with a number 
of previous authors we suggest that the inner regions of the accretion discs in FU Ori 
objects are evacuated by the action of a magnetic propeller anchored on the central 
star. We develop an analytic solution for the steady state structure of an accretion 
disc in the presence of a central magnetic torque, and present numerical calculations 
to follow its time evolution. These calculations confirm that a recurrence time that 
is consistent with observations can be obtained by selecting appropriate values for 
viscosity and magnetic field strength. 

Key words: accretion, accretion discs - stars: FU Orionis - stars: T Tauri - stars: 
magnetic fields - stars: pre-main-sequence 



1 INTRODUCTION 

; The young stellar objects (YSOs) are stars which have 
' not yet completed the process of star formation. They fre- 
quently retain some portion of their proto-stellar discs, 
from which mass accretion is thought to take place. The 
T Tauri stars are YSOs with masses < 2 Mg, and 
accre tion rates in th e regio n ~ 10~ 7 — 10 -8 Moyr -1 , 
' e.g. lHartmann et al.1 <ll998f) . T Tau stars are divided 
into two subclasses, Classical T Tau stars (CTTs) and 
Weak Line T Taus (WTTs). CTTs have been shown to 
posses discs by the observation of double peaked absorp- 
tion lines, e.g. lHartmann fc Ke nycvn|j[l^ ^) . and from the 
presence of an infra-red ex cess ijAdams. Lada fc Shulll987t 
iKenvon fe Hartmannlll987l) . CTTs have also been observed 
to produce extended emission flBeckwith et al.lll984T) which 
is consistent with th e existence of an acc retion disc. WTTs 
are a faster spinning feouvier et al.ll995l) type of YSO with 
weak absorption lines and a less pronounced infra-red excess 
( Hart mannlll998h . The presence of a disc is less certain in 
these systems. 

Some YSOs have been observed to undergo outbursts, 
during which they increase in brightness by up to two or- 
ders of magnitude on a timescale t T i sc ~ 1 yr. These ob- 
jects, termed FU Orionis stars, are otherwise very simi- 
lar to CTTs. Therefore it is generally assumed that FU 
Ori stars are members of the same population as T Tau 



* email address: owen.matthews@astro.le.ac.uk 



stars, e.g. lHartmann fc Kenvonl lll985l) . The outburst du - 
ration has been estimated as t b ~50 yr ijHartmannl f 1 998 \) . 
from the decay times of the outburst light curves (e-folding 
decay times of the order of decades are typical). Both 
t b and the recurrence time of the outbursts (t rcc ) re- 
main poorly known as no stars have yet been observed 
to return to their pre-outburst states. The peak outburst 
accretion rate is ~ 10 -4 Mq yr -1 , which is found by 
matching numerical simulations to outburst spectra, e.g. 
IKenvon. Hartmann fc Hewettl jl988l) . The mean accretion 
rate during outburst is ~ 1CP 5 Mq yr -1 . For a detailed re- 
view of the FU Ori and T Tau stars, and re ferences, see 
iHartmannl lll998f) : lHartmann fc Kenvonl il996h . 

The outbursts of YSOs are not yet fully understood, 
but are believed to be associated with a thermal-viscous ac- 
cretion disc instability. Models based on the thermal-viscous 
disc instability have widely been used as an explanation for 
similar outbursts in other stell ar systems, esp ecially binary 
stars such as dwarf novae, e.g. IWarnerl |)199!tI) . However the 
long recurrence times of YSOs do not correspond with those 
predicted by the standard thermal-viscous disc instability 
model (DIM) with an accretion disc viscosity of the same 
order as that applied to binaries. In this paper we argue that 
the standard DIM, with a quiescent Shakura-Sunayaev vis- 
cosity parameter of a c > 10 -3 , in the case of an accretion 
disc with an inner radius comparable to the radius of the 
central star would predict a series of shorter, more frequent 
outbursts than are observed in any YSOs. One possibility is 
that the observed outburst behaviour could be reproduced 
simply by decreasing the Shakura-Sunyaev viscosity param- 
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eter to a very low value e.g. iBell fe LirJ Jl994f) . We propose 
an alternative model which incorporates the effect of a stel- 
lar magnetic field on the inner accretion disc. A magnetic 
field anchored on a rapidly rotating star can act as a pro- 
peller (e.g. Lovelace, Romanova & Bisnovatyi-Kogan 1999; 
Wynn, King & Home 1997) creating a low density region 
in the inner disc, and reducing accretion to the order of the 
CTT rate. Outbursts are thought to begin when the sur- 
face density exceeds a critical value E cr §. anywhere in the 
disc. Sinc e the critical density increases with distance from 
the star iCannizzo. Shaffer fc Wheeler! fl988l) . a truncated 
accretion disc is able to accommodate more mass, forcing 
outbursts to be more dramatic and much less frequent. 

In Section |5| we review estimates of the recurrence and 
rise times of outbursts in T Tau systems. These estimates 
are compared with those expected from the standard DIM, 
we also outline the effect of viscosity. An extension of the 
DIM, including a magnetic propeller, is discussed in Section 
An expression for the evolution of the disc in the radial 
dimension in the presence of a magnetic field is derived in 
Section|3]and an analytic solution is obtained for the steady 
state. Finally numerical calculations are performed which 
demonstrate the feasibility of using a magnetic propeller to 
create the required depleted region in the centre of a T Tau 
disc. 



2 OUTBURST TIMESCALES 
2.1 Disc instability model 

The thermal- viscous disc instability model is now generally 
accepted as the mechanism behind outbursts in, for example, 
dwarf novae, e.g. (Buat-Menard, Hameury & Lasota, 2001; 
Warner, 1995), and other accreting binary systems. However 
The standard DIM is unable to explain the long timescales 
associated with the outburst cycle of FU Ori stars at least 
for similar viscosities to those applied to discs in binaries. 

In the thermal- viscous instability model a limit-cycle in 
temperature and density forces the accretion disc to cycle 
between two stable states: a hot, high accretion, high vis- 
cosity state; and a cold, low accretion, low viscosity state. 
Once an annulus of a disc reaches the hot state it triggers a 
heating wave leading to an outburst which affects the whole 
disc, or at least a portion of the disc which extends to some 
maximum radius. Th e DIM is discussed in detail by e.g. 
iFrank. King fc Rainel l)2002h : ICannizzo. Shaffer fc Wheeled 
( 1988). The trigger causing the change from the cold to the 
hot state can be expressed in ter ms of a critical value of the 
disc surface density E. Following ICannizzo et all iTlfi 
approximate this critical E with the formula 



E|?;f t h = n.4^0 05 Mr°- 3 V 



-0.86 -2 

g cm 



(1) 

where Rio is the radius from the star in units of 10 10 cm, Mi 
is the mass of the star in solar masses, and a c is the Shakura- 
Sunyaev alpha viscosity during the cold (quiescent) state. It 
is clear that the lowest E c ** , which would trigger outbursts, 
must appear close to the inner edge of the disc (Ri) and that 
outbursts would begin near Ri unless some process were to 
cause a distortion of the disc's surface density profile. The 
radius at which the outburst is initiated has a strong effect 
on both t D b and tree- 



It is possible to estimate the outburst duration t b by 
considering the e-folding decay time t c of the observed light 
curves of FU Ori stars. These timescales are typically in 
the range of lOyr < t e < 50yr. An FU Ori star in outburst 
is about two m agnitudes brighter than in quiescence, e.g. 
iHartmaiml l)l998h . If the outburst is considered to be over 
when the star returns to its quiescent luminosity then as a 
first approximation we may say that the outburst time is of 
order t b ~ 0.4i e Amln 10 <10 2 yr, where Am is the change 
in magnitude between outburst and quiescence. 

Since no FU Ori stars have been observed through a 
complete outburst cycle it is difficult to estimate t ICC from 
observational data. However it is possible to obtain a rough 
estimate of t rec by comparing mean mass transfer rates in 
quiescence, (M qu ), and during outburst, (M b). Assuming 
that a large proportion of the region of the disc involved in 
the outburst is accreted, this mass must be replaced in a 
time ircc and hence 

{Mob)t oh ~ (M qu )t re c • (2) 

Solving relation J2J yields t TCC ~ 5000 yr, where we have 
adop ted (M au ) ~ 10"^ Moyr -1 and (M ob ) ~ 10 -5 Meyr -1 
from IHartmannl ljl998l) and the estimate t h ~ 50 yr. This 
estimate is likely to be an upper limit as the mass transfer 
rate through the disc in quiescence must be larger than the 
observed accretion rate to allow an accumulation of mass 
leading to the next outburst. 

A further approximate analytic limit for t lac in the DIM 
can be obtained by considering the mass flux through the 
annulus in the disc at which the outburst is initiated. The 
surface density of this region must reach E c ^. in a time t rcc , 
leading to the relation 

M(Ri) 271-fliEA.R 



M(Ri) 



(M qu ) 



(3) 



where M(Ri) is the mass contained in the annulus, which be- 
gins close to the inner radius Ri of the disc and extends out- 
wards a distance AR. Approximating AR ~ 0.1 R (which is 
not unreasonable if an aspect ratio of H/R ~ 0.1 is adopted 
throughout the disc) it is possible to express t rcc in the form, 



(M qu > 



(4) 



During quiescence the surface density in the disc is limited 
by E < E^f t h . Therefore equations © and JTJ can be com- 
bined to give 

yr , (5) 



tree < 2 x 10- 2 Mf - 36 _Rf; 1 ° 5 (M qu , 15 )- 1 Q ! - 1 



where R u io represents the inner radius of the disc in units 
of 10 10 cm, (M qu , 15} is the mean quiescent accretion rate 
in units of 10 15 gs _1 . In this paper we adopt a value of 
a c ~ 0.01 which is the quiescent viscosity parameter com- 
monly applied to dwarf novae accretion discs. This assump- 
tion will be discussed further in section 12.21 The estimate 
of the surface density profile that is used here is of course 
an upper limit, it will be seen however in Section 0] that 
this is a reasonable approximation. The estimate © can 
be tested by comparing its predictions to the outburst 
behaviour of the dwarf novae, which have been observed 
through many outburst cycles. In a dwarf nova the accret- 
ing star is a white dwarf with a mass M+ ~ 1 Mq, and 
radius i?* ~ 0.05 Rq, and the quiescent accretion rate is 
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(AV} < 1CT 11 M yr _1 (IWarnerlll995l) . Using these values, 
and assuming that Ri ~ i?*, equation gives an estimate 
for the recurrence time of t rcc ~ 10 d, very close to the 
observed t ICC in these systems. In order to apply to a 
typical FU Ori system we adopt values of M* ~ 1 Mq, and 
(M qu ) ~ 1(T 7 M yr _1 . If the disc is assumed to extend 
down to a typical CTT stellar radius of Ri ~ 2Rq, then 
equation JSJ yields a recurrence time of i rec ~ 0.5 yr This 
is clearly far too short to be consistent with observation. 
Indeed, to obtain t rcc ~ 5000 yr we require an inner disc 
radius of Ri ~ 40 Rq. 

The timescale on which the luminosity of the disc in- 
creases (irise) in an FU Ori outburst is governed by the ther- 
mal timescale of the region in which the outburst begins , 
in the case of an outside-in outburst. Clar ke fc Sverl (1996) 
have examined this is detail and point out that rise times of 
<1 yr imply that outbursts must begin at radii much larger 
than that of the central star. They estimate values >15 Rq, 
in rough agreement with the estimates above. If we assume 
the observed values of t r i se are of the same order as the local 
thermal timescale ith at th e point where the o utburst begins 
t t h we obtain the relation l|Frank et al.ll20o3 l 



irisc ~ t th (R) ~ (an k (R)) 1 



a c \GM* 



(6) 



where flfc (R) represents the Keplerian angular velocity at 
a given radius. The alpha viscosity parameter in quiescence 
is once again assumed to be a ~ 0.01. The observed rise 
time of £ri sc ~ 1 yr is consistent with an inner disc radius 
of -Ri ~ 40 R©, which is similar to the estimate obtained 
from t r ec and consis tent with the large radii predicted by 
IClarke fc SveiUl996l) . 



2.2 Low viscosity models 

It has been shown that the above estimates for t Tec and irise 
are consistent with the DIM if it is modified, for example, 
by including a truncated inner disc where Ri ~ 40 Rq. We 
note here however that long values of t rcc and t b, consis- 
tent with observation, can be reconciled with the standard 
DIM if a low value of the quiescent disc viscosity is invoked. 
Equation © has a dependence on a c and so reducing a c 
to small values can yield a t lac which i s arbitrarily long . 
Ind eed such a model has been applied bv lBell fc Linl l|l994h 
and lBell et al.l IJ1995T) who use a one dimensional disc m odel 
similar to that described in SectionH iBell fc Linl <ll994l) sim- 
ulate full outbursts and include detailed thermodynamics in 
order to calculate the onset of the thermal instability. This 
is then associated with an arbitrary increase in the disc vis- 
cosity. For most of these simulations E = at the inner 
boundary and M is held constan t at the outer boundary. 
With this model IBell fc Linl Jl994T) are able to obtain recur- 
rence times of t ICC ~ 1000 yr using the standard DIM with 
a c ~ 10 -4 . Our analytic estimate of t Tec from Equation 10 
requires an even lower value of a c ~ 10 -6 to obtain t TCC of 
this order. This discrepancy is because the analytic estimate 
developed in 12. II includes the assumption that the outburst 
begins at a radius of R = Rj ~ 2Rq whereas according 



to the resuhs_c4jBell&^in| l)l994T) the outburst begins at 
R ~ IORq . IBell fc Linl (|l994l) also include a feedback mech- 
anism whereby the radius of the disc is a function of M. 



We note here that this model is successful, providing that 
outbursts can be initiated at -R > 10 Rq. Of course a simi- 
larly long t rcc could be replicated by outbursts initiated at 
smaller radii assu ming even lower viscosities. Nonetheless 
IBell fc Linl ( 1994) do confirm that the dwarf nova like es- 
timate of a c ~ 10~ 2 is far too high to reproduce values of 
t ICC in agreement with observational limits in the standard 
DIM. 

Theoretical work has been done which suggests that 
a very low value of a c is possible in protostellar discs, 
because the region in which thermal ionisation is able to 
occur throug hout the thick ness of the disc is limited to 
< 0.1 AU jGammielll996T) . However the FU Orionis out- 
bursts are driven in the inner disc, which is likely to be ther- 
mally ionised, and where as a resul t the viscosity is likely to 
be similar to that in dwarf novae. lHartmann et al.l (ll998T> 
have argued that the radii of circumstellar discs can be used 
to constrain the quiescent viscosity parameter to a c ~ 0.01, 
similar value to the dwarf nova case, throughout the disc. 
The remainder of this paper will be dedicated to the in- 
vestigation of an alternative mechanism which can extend 
recurrence times, by disc truncation, to within the observa- 
tional limits while using a dwarf nova like low state viscosity 
of a c ~ 0.01 in the inner disc. 



3 DISC TRUNCATION BY A MAGNETIC 
FIELD 

In Section 2 we argued that the standard thermal-viscous 
instability outburst model is insufficient to explain the out- 
bursts of the FU Ori stars. A truncated disc however pro- 
duces reasonable estimates of t lac and i r i sc . One possible 
mechanism for truncating the disc is a magnetic field asso- 
ciated with the accreting star. 

The suggestion that the accretion disc in YSOs is trun- 
cated by the magnetic field of the young star, or other- 
wise, is not new and has appeared in the literatur e man y 
times. Examples include lUltchin. Wvnn fc Reeevl feOOll) . 
who model the CTT infra -red excess with a truncated disc, 
lArmitage fc Clarke! l)l99ffl who expl ain T Tau spin regula - 
tion with a truncated disc model, and lClarke fc Sverl f 1996). 
A magnetic mechanism which may lead to disc truncation 
is investigated below. 



3.1 The magnetic propeller 

The magnetic fields of the T Tau and FU Ori stars are 
widely believed to play an important role i n their evolu - 
tion and behaviour, e.g. lYi fc Renvoi] (ll997l) : ISafieil <199Sft . 
Good estimates of the magnitudes of these fields are, how- 
ever, difficult to obtain. The most recent estimates available 
of the magnetic field of t he CTT star BP Tau, measured 
using spectropolarimetry d.Tohns-Krull et al~lll999F) . put the 
value at B — 2.4 ±0.1 kG. This is in good agreement with 
earlier measurements for the magnetic fields of both CTTs 
and FU Ori stars, which ha ve been calculated using Zee - 
man broadening effects, e.g. I.Iohnstone fc Periston! l)l987l) . 
The inner accretion disc in a CTT will be partially ionised, 
as it is irradiated by the accreting star, and will therefore 
interact with the stellar magnetic field. Angular momentum 
can be transfered from the disc to the star and vice-versa, 
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via this magnetic interaction. The magnetic field is therefore 
important in determini ng the spin rates of these stars as dis- 
cussed in more detail in lPophaml (|l99(J l: lArmitaee fc Clarke) 
(1996). 

The magnetic propeller effect (e.g. Wynn et al. 1997; 
Lovelace et al. 1999) involves a net transfer of angular mo- 
mentum from the star to the disc via the magnetic field. In 
this case the interaction between the disc and the field takes 
place at radii larger than the co-rotation radius R co , the 
point at which the circular Keplerian angular frequency Qk 
of the disc matches the angular frequency fl* of the star. For 
the purposes of this paper we assume that the stellar field is 
dipolar and that it rotates at the same rate as the star such 
that 



Ok(R c 



(7) 



where Qb represents the angular frequency of the magnetic 
field. At R co it is clear that a body in a circular Keplerian 
orbit remains stationary in the frame of reference of the 
magnetic field, and there is no torque between the disc and 
the field at this point. Hence, the magnetic timescale t mag 
(the timescale on which magnetic torques redistribute an- 
gular momentum in the disc) is infinite at R co . Outside R co 
we have £Ik < fi* and the disc gains angular momentum at 
the expense of the star. Matter which is accelerated in this 
manner will be pushed out to greater radii until the mag- 
netic torque (which for a dipole field decreases as ~ R~ 6 ) 
becomes small when compared to viscous stresses in the disc. 



This occurs when t n 



> U 



where t v - lsc is the timescale 



on which viscous stresses redistribute angular momentum 



within the disc. The radius at which t n 



U 



is usually 



termed the magnetic or magnetospheric radius i? mag . Inside 
Rco, the magnetic torque extracts angular momentum from 
the disc. In this case the disc mass is rapidly accreted onto 
the star. 

The results of both the propeller and the accretion 
regimes are the formation of a depleted region in the disc 
close to the star. The radius of this region is Ri ~ i? mag , 
and therefore is a function of B, the spin period P sp i n of 
the star, and the disc viscosity. If we assume, as argued in 
Section |21 that an inner disc radius Ri ~ 40 is required 
to be consistent with observational constraints then Ri must 
exceed R co even for slow spinning T Taus. Although both 
the accretion and propeller regimes lead to the formation of 
a depleted region in the centre of the disc, only the propeller 
mechanism leads to the accumulation of disc mass at radii 
greater than R co - 

Below we derive an expression for the magnetic 
timescale t mag in a completely magnetised accretion disc in 
the presence of a dipolar ma gnetic field. The result is almost 
identical to that adopted bv lLivio fc Pr inglc ( 1992). The in- 
clusion of a magnetic field in a hydrodynamic system in- 
troduces two additional terms to the Euler equation; a mag- 
netic p ressure term and a magnetic tension term, e.g. iDendvl 
( 1990). The magnetic pressure term is negligible where B is 
small. The magnetic tension term can be expressed as 



1 

PAVc 



(8) 



where p<j represents the density of the disc, r c is the local 
radius of curvature of the field lines and B z and B v repre- 
sent the vertical and azimuthal components of the magnetic 



field respectively. We use the a pproximations r c ~* H B z j 
JPearson. Wvnn fc Khi3ll997l) . and H ~ O.li?. The ratio of 
field strengths can be expressed in the form (e.g. Livio & 
Pringle 1992) 



Btf 

~bI 



(n k - 



(9) 



The magnetic timescale can be defined in terms of the mag- 
netic acceleration and the Kepler velocity by the relation 



tn 



2nR 2 p d 



n't 



5Bi (f2 k - Q*) 



(10) 



Finally the volume density can be related to the surface 
density by the expression pa ~ T,/H ~ lOE/i? and, for a 
dipole field, we have B z ~ fj,R~ 3 where p is the magnetic 
moment of the star. These relations lead to the following 
expression for the magnetic timescale 



tn 



47rE % /GAO? 11/2 



([&: 



3/2 



(11) 



It is of course possible to use other prescriptions for the mag- 
neti c timescale, for exam ple the diamagnetic case as used by 
e.g. IWvnn et al. I il99?fl . In general however £ mag can nor- 
mally be parametrized in the form 



tn 



2Y, 



P 



{[■&: 



3/2 



(12) 



where 7 and /3 are parameters determined by the magnetic 
interaction model. We use this general expression to derive 
equation 12511 in Section 0] In the fully magnetised case rep- 
resented by equation <f 1 1 pi we have 7 = 7/2 and is deter- 
mined by the mass and magnetic moment of the star accord- 
ing to 



(13) 



3.2 Magnetic models for CTTs and WTTs 

Several authors have co nsidered magnetic truncation of the 
accretion disc in YSOs. iKenvon. Insu fc Hartmannl (|l996h 
use this idea to model the infrared excesses of CTTs. 
lArmitaee fc Clarke! l)l99fiF) model a magnetically truncated 
accretion disc and use this model to explain the regulation 
of the spin periods of T Tau stars on a timescale ~ 10 yr. 
The magnetic propeller effect has pre viously been applied t o 
the case of T Tau accretion discs bv llJltchin et al.l l)200ll) . 
who use a smoothed particle hydrodynamics technique to 
show that a central hole is formed by the magnetic field, 
exactly as described above and as will be confirmed by nu- 
merical results in Section^] In the full three dimensional case 
some accretion should occur, as there will be a component 
of the orbital velocity which is not in the azimuthal direc- 
tion with respect to the dipole field. The accretion stream 
follo ws the magnetic fi eld lines of the accreting star. Indeed 
ICullbring et al.l l)l998t) find that CTT spectra agree with an 
accretion stream of this sort better than they do with stan- 
dard boundary layer accretion. This is analogous, in binary 
systems, to the flow in the case of t he intermediate pola r 
(IP) stars as discussed for example in IMurrav et al.Ul99Sft . 
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iLivio fc Pringld il992l) discuss the effect of a magnetic pro- 
pelier on the onset of outbursts in dwarf novae. 

Magnetic models have also been used to explain the 
transition between CTT and WTT states by allowing the 
magnetospheric radius -R mag to move, which affects the ac- 
cretion behaviour of the system. When i? mag is large, the 
disc is propelled outwards to large radii and very little ac- 
cretion takes place. This state can be identified with a WTT 
system since there would be a low accretion rate and no in- 
ner disc. If the star is in a state where 7? mag is smaller, then 
it may be that _R mag < -Rco, in which case the system would 
become a magnetic accretor. As an alternative to the mag- 
netic accretor, a weak propeller could be invoked. In this case 
Rmag > Rco but -R ma g is still smaller than the strong field 
value. In either case these systems could then be identified 
as CTTs as a disc would exist close to the star, although 
not extending to the surface. This disc could produce the 
observed infra-red excess and emission line characteristics 
of CTTs. 

IClarke et al.l 1^995) suggest a 'magnetic gate' where 
magnetic cycles in the star cause a change in the magneto- 
sp heric radius, R mag . A secon d possibility is that proposed 
bv lArmitaee fc Clarke! 1^)96!) where i? mag can be changed 
by spinning the star up or down by disc interactions. High 
spin stars, with a greater i? mag , would then be classified as 
WTTs, which is in agreement with observation. A rapid spin 
up can most simply be achieved by a high accretion phase 
such as an FU Ori outburst. The gradual spin down of the 
system in its W TT phase and trans i tion t o CTT is mod- 
elled in detail bv lArmitage fc Clarke! ^)<Mi). The effects of 
a magnetically truncated disc on the spin evolution of the 
T Tau stars will be considered in Section 13.41 



3.3 A magnetic model for FU Orionis outbursts 

As discussed in Section 13.21 the action of a magnetic pro- 
peller in the disc of a T Tau object will cause the truncation 
of the disc at some point outside R co . The precise position of 
truncation depends on the strength of the stellar magnetic 
field, the rate of spin of the star and the viscosity of the disc. 
If a depleted region, or hole, is formed then the total mass 
Mob of the disc is allowed to reach a higher value before 
E'"? t h (RC) is exceeded and an outburst begins, increasing 
both t 00 and t rcc . During the quiescent, or propelling, phase 
mass transfer is not completely prevented, but is reduced 
to the CTT accretion rate of M ~ 10" 7 M Q yr _1 . Some ac- 
cretion will still occur because in three dimensions, as men- 
tioned in Section f3. 21 some component of the azimuthal ve- 
locity is not perpendicular to the magnetic field. The mass 
which leaks through the magnetic barrier accretes along the 
magnetic field lines onto the central star, which is analo- 
gous to the behaviour of other magnetic propeller systems. 
In a forthcoming paper we attempt to apply a similar mag- 
netic mechanism to the cataclysmic variable WZ Sagittae, 
in which case the quiescent accretion rate is ~ 1 percent of 
the outb urst rate. Indeed ev en in the interacting binary AE 
Aquarii (|Wvnn et al.lll99/ft . a propeller system so strong 
that it completely prevents the formation of a disc, the ac- 
cretion rate is ~ 1 percent of secondary mass loss rate. It 
should be noted that this accretion is an intrinsically three 
dimensional effect and that the situation is somewhat differ- 



ent in the one dimensional model applied later in the paper. 
This distinction is discussed in section |4~T1 

It is of course possible to imagine a disc in which the 
propeller is sufficiently strong to prevent outbursts. In this 
case a steady state may exist in which the accretion rate 
is constant and the critical density is never reached. The 
parameters for which this is possible will be investigated 
in section 14.21 It should be stressed here however that in 
order for any steady state to exist, some 'leakage' of mate- 
rial through the magnetic barrier is required. The relative 
frequency of outbursting and non-outbursting CTTs will de- 
pend on the distribution of magnetic moments and of disc 
viscosities. The discussion of this effect requires a full pa- 
rameter search and is left to future work. 

If the disc goes into outburst, the viscosity is believed 
to increase by up to an order of magnitude. The resulting 
decrease in i v isc reduces R ma g to a point much closer to, or 
even within, the surface of the star, and causes M to increase 
by several orders of magnitude to the peak FU Ori accretion 
rate of M ~ 10 -4 MQyr -1 . The dissipative (viscous) effect 
dominates over the advective (magnetic) effect at all radii 
in this case. This continues to be the case until much of the 
stored mass in the disc has been accreted and the disc is 
allowed to return to the cold state. Indeed for as long as 
the condition Ri < R co is satisfied, the magnetic field will 
actually assist accretion in the inner disc. 

On the disc returning to the cold state, Ri is forced out 
once more to its original position outside R co , assuming that 
neither B nor P ap in are significantly altered by the outburst. 
The system should now return to its pre-outburst state, with 
a CTT accretion rate, and the mass of the outbursting region 
should begin to build up once more. Hence, providing that 
a magnetic propeller can produce a depleted region, it can 
be invoked to explain the long t lac as well as the long J r i se 
in FU Ori systems. 



3.4 Spin evolution 

Since P S pm effects R mag and hence has an influence on t lec , 
it is important to consider the effect of spin evolution on 
FU Ori outbursts. This can be done to a first approximation 
by comparing the spin up timescale of the star to t rec and 

tab- 

According to lArmitaee fc Clarkd jl996t) the squared ra- 
dius of gyration of a T Tau star is k 2 ~ 0.2 for most of its 
lifetime. We adopt this value of k 2 and assume a solar mass 
star with a typical T Tau spin period of P sp m ~ 3 d. The 
rotational angular momentum of the star is given by 

Li, = ^M+Rttt* . (14) 

Away from the inner radius angular momentum transport in 
the accretion disc is dominated by viscous effects. It is only 
at the inner edge of the disc that magnetic effects become 
important. 

During quiescence the magnetic propeller will cause 
the star to experienc e a net spin down torque. 
lArmitage fc Clarkd dl996ll suggest that a propeller will spin 
down the star on a timescale of i sp in,qu ~ 10 5 yr. It is likely 
then that stellar spin evolution is insignificant during a sin- 
gle quiescent phase. 

In the case of magnetic accretion the spin up torque on 
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the star is given by = Ml(Ri), where I (Ri) represents 
the specific angular momentum of matter at Ri, giving 

Li, = My/GNLR]' 2 . (15) 

The spin up timescale can be expressed as 

i spin ~ L t /U . (16) 

During outburst (M oh ) ~ 1 X 10" 5 M©yr _1 and Ri ~ i?* ~ 
2 R0 because the central region is likely to have become 
filled. From equation 1)16)1 . this yields t sp i n ,ob ~ 10 3 yr. 
Therefore even in outburst £ S pin,ob > *ob compared with an 
estimated t a b ~ 10 2 yr. In this case it is sensible to assume 
that the measured M is the same as that at the surface of 
the star and to neglect transport of angular momentum from 
the star to the disc. This result shows that spin evolution 
should therefore not have a significant influence on outburst 
behaviour during one outburst, although there may be a 
measurable cumulative effect. 

Although spin evolution may turn out to be significant 
in the longer term for T Tau stars, it is still possible to 
illustrate the principle of outburst and recurrence time be- 
haviour with a fixed spin and hence a fixed -Rmag- It is as- 
sumed for the purposes of this paper that the outbursts oc- 
cur in the CTT phase and that the system returns directly to 
a similar CTT state after outburst. Any magnetic evolution 
of the T Tau star is neglected. 



4 THE STRUCTURE OF A TRUNCATED DISC 
4.1 Analytic Solution 

A thin, axisymmetric accretion disc is assumed to be ro- 
tating about a central mass M*. As we are only concerned 
with the inner Sections of the disc where outbursts begin, 
and the total dis c mass can be estimated to be ~10 -3 M© 
iflartmannlll99st) . it is reasonable to neglect self-gravity in 
this portion of the disc. Averaging the basic hydrodynamic 
equations over all azimuthal phases and integrating over the 
direction normal to the plane of the disc, the continuity 
equation becomes 

gE I BjBBvB) 

dt R dR ' 1 ' 

the radial component of the Navier Stokes equation becomes 



dvR 
dt 



+ Vr 



d 



'VR 

OR 
d 



3i? 3 / 2 dR 



R 



R 3/2 vT, 



dp 
dR 



dvs, 
dR 



,GM 

d{R 2 vT,VR) 



(18) 



3R 3 



dR 



and the azimuthal component of the Navier Stokes equation 
becomes 



dv v vr d(Rv v ) 



dt 



+ 



R dR 



R 2 dR 



(19) 



where v v and vr denote the azimuthal and radial compo- 
nents of the velocity respectively and v represents the kine- 
matic viscosity of the disc. 

By multiplying the azimuthal component 11911 of the 
equation of motion with R, we obtain an equation for the 
specific angular momentum in the disc, 



'dl dl\ 
M +VR dRj 



1JL 

RdR 



i? 3 i/E 



_d_ 

dR 



+ EA 



(20) 



Here, we have already assumed an additional external torque 
acting on the disc, where A is the injection rate of angular 
momentum per unit mas s, following the same procedure as 
iLin fc Papaloizor |l986). A similar procedure is also em- 
ploved by lPrhigk (m)9lT) in the case of a binary system. 

Assuming a sufficiently cold disc and therefore that the 
dynamical timescale tdyn ~ R/v v is much smaller than any 
other timescale, we can adopt a Keplerian approximation 
for the azimuthal motion such that v v = (GM/ R) 1 ^ 2 . Then, 
12m can be solved for the radial velocity in the disc, 



Vr = 



0(*Va„E) | 2A R^ 2 



dR 



(21) 



Inserting equation 12 It into the continuity equation 1)170 re- 
sults in an evolution equation for the surface density of the 
disc, 



<9E _ 3__d_ 

dt RdR 



R 



1/2 



d{R 1/2 vT.) 
dR 



1JL 

RdR 



2AE 



r3 /2 



,(22) 



where the right hand side is composed of a diffusion term 
on the left and an advection, or torque term on the right. 
To parametrize the specific torque A, we write 

th t\ 

where t\ is the timescale on which the local disc material 
gains angular momentum. 

In the present case, the source of the torque is the mag- 
netic interaction of a rapidly rotating, moderately magnetic 
YSO within the disc. The torque timescale is therefore equiv- 
alent to the magnetic timescale, th ~ imag- This magnetic 
timescale is given in equation 112H . 

In general, the kinematic viscosity v may depend on 
radius as well as on surface density, v = v(R, E). For the 
purposes of this paper v is taken to be independent of E and 
parametrized according to a power law such that v = 5R e . 
Kinematic viscosity may also change with time according 
to an outburst limit cycle. It may be possible to implement 
this form of v and, with suitable boundary conditions, to 
simulate complete out burst cycles in the future. The viscous 
timescale is given by jFrank et al. ^ 2002) 

By considering equation 112H . a general expression for the 
advection term in equation is obtained 



I JL 

RdR 



2AE 



R 3 / 2 



RdR 



1 



R 



Ro 



3/2 



(25) 



Equation 122H can be solved analytically for the steady state. 
By definition in the steady state M = — 27ri?EuR, is a con- 
stant throughout the disc. Using equation 1)21)1 and with a 
magnetic torque as defined in equation 1231 this condition 
can be used to derive a general solution for E in steady state, 



3/2 (2-7) 



(26) 



_ M / r_R 

3nS + 3(2- 7 )5 [Ir co 



where C is an arbitrary constant. As will be seen this equa- 
tion has a form such that E becomes negative close to the 
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star. The inner part of this steady state solution is thus 
clearly unphysical. It is difficult therefore for a physically 
motivated inner boundary condition to be chosen. Further- 
more the general shape of the function is rather sensitive to 
this boundary condition. The outflow condition dT,/dR = 
will be applied however as a first approximation. The inner 
boundary is taken to be very close to the surface of the star. 
This gives, in its most general form, 



C = 



PR*~' ( (7 + ff i/2 , (I-7-e) Rj \ . . 



eMR 



1/2 



37r(| + e) ■ 

This analytic solution for the steady state is precisely repro- 
duced by the numerical solution in the following Section. It 
is however instructive to examine the analytic form of the 
solution. It can be seen from the form of equation 1261 that, 
counter- intuitively, in steady state the position of Ri, defined 
here as the point where E = 0, is independent of S. The disc 
structure at very large radii tends to the non-magnetic case 
as would be expected. The non-magnetic solution is depen- 
dent only on M, and v. It is almost ide ntical to the solutio n 
for the steady thin disc presented bv iFrank et al. I {2002). 
Indeed this solution, 



M 
3tt 



I - 



R 



1/2 



(28) 



can be reproduced precisely by applying an inner boundary 
condition such that E = 0. In this case, as well as in the full 
magnetic case it holds that E oc 1/u. 

The steady state solution in one dimension cannot rep- 
resent a physical situation. In steady state M is constant 
throughout the disc by definition. However in the steady 
state solution found above there is an evacuated region in 
the disc through which no matter can pass. Numerically this 
represents the transport of negative mass, but this clearly 
has no physical meaning. In three dimensions however it is 
possible to imagine similar behaviour by allowing mass to 
flow along magnetic field lines in the inner disc. The one di- 
mensional solution can then be taken as an approximation 
of the behaviour of matter in the plane of the disc. This 
approach is given some validity by the similarity of these 
one dimensional results to those obtained i n three dimen- 
sional simulations. This can be seen in e.g. lUltchin et al. I 
(2001) and also in our own preliminary three dimensional 
computations. 



4.2 Numerical Solutions 

It will be shown by solving equation (1221 numerically, that 
it is feasible to create a depleted region of Ri ~ 40 Rq in the 
centre of a T Tau accretion disc with a magnetic propeller 
and hence to increase t TCC to arbitrarily long times. The 
code is similar to that used by lArmitaee fc Clarke! dl996l) 
to model a T Tau disc about a magnetic star, and their re- 
sults are similar to those p resented below . The results are 
also comparable to those of IPringlel dl99ll) . In all cases the 
simulations tend towards a steady state which can be ex- 
pressed by the analytical solution derived above. The ana- 
lytic solution has been deliberately kept in as general a form 



as possible, howeve r for the numerical calculations e = 3/4 
(Fran k et al.ll2002T) and 7 = 7/2 have been adopted. The 
computation is sensitive to resolution. Convergence with the 
analytic solution only occurs at high resolution and when the 
inner boundary conditions match. 

The solver code evaluates the diffusive and advective 
terms separately and combines them using operator split- 
ting. The diffusive term is computed with a modified Crank- 
Nicholson scheme jPress et al.lll992h and the resulting val- 
ues of E are use d as input for a m odified two-step Lax- 
Wendroff scheme ijPress et al.lll992h which is used to solve 
the advective term. 

The first pair of simulations were performed starting 
from a Gaussian density distribution, using the code de- 
scribed above. This initial distribution represents a viscously 
spreading ring. One simulation was with, and one without, 
an advection term, representing a magnetic field. The pa- 
rameters used in the code were as follows: the mass of the 
accreting star was taken to be M* = 1 Mq . The co-rotation 
radius was calculated on the basis of a star spinning with a 
typical T Tau period of 3 days which yields R co ~ 10 Rq. 
The boundary conditions applied at both the inner and outer 
radii of the disc are outflow conditions and the initial den- 
sity distribution was centred at a radius of 5 -R co . The disc 
was simulated from 0.1 R co ~ 1 R© to 10 R co - The value 
(3 ~ 10~ 13 is consistent with an object of R ~ 1 R© and 
B ~ 1 kG. 

With the parameters described above, the effect of an 
advective term on the surface density is illustrated in Fig- 
ure In the non-magnetic case, where the diffusion term 
acts alone, the initial Gaussian mass distribution quickly 
moves towards the accreting star, as expected. Thi s is very 
similar to the well known result of IPringlel il98ll) . In the 
magnetic case, the evolution of the outer edge of the mod- 
elled region is very similar to the non-magnetic case. This 
is not surprising considering the R~ 7 ^ 2 dependence of the 
advection term and is in agreement with what was expected 
from the analytic solution. At the inner edge however, it can 
clearly be seen that a depleted region of radius ~ 1.5 R co 
is formed. It is interesting to observe how matter in the 
initial Gaussian distribution is transported by the diffusion 
and advection terms. Radial mass transfer rate M is plotted 
in Figure The magnetic and non-magnetic cases appear 
quite similar upon first inspection. Initially the movement 
of the mass in the Gaussian is dominated by a diffusion ef- 
fect. This is represented by the two large lobes. However, in 
the non-magnetic case the flux near the inner boundary in- 
creases in magnitude with time. This represents an accretion 
process. In the magnetic case there is no accretion, and the 
inward flux rapidly disappears with time and the outward 
flux, driven by advection, is dominant at later times. 

Further simulations were performed using a constant 
inflow of M ~ 10~ 7 M©yr _1 at the outer boundary and 
an outflow condition at the inner boundary. In this case 
an empty disc is assumed at the beginning of the simu- 
lation. These conditions are qualitatively a more realistic 
representation of the formation of the inner disc after an 
outburst, and should more correctly illustrate where mass 
build up would first trigger a new outburst. Figure shows 
the outcome of these simulations with varying magnetic 
field strengths. The critical density Ej."? t h is marked on these 
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(a) 



0.02 Viscous Times 
0.10 Viscous Times 
2.00 Viscous Times 



(b) 



S 3e-05 - 
< 



0.02 Viscous Times 
0. 10 Viscous Times 
2.00 Viscous Times 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 
R7Rc„ 



Figure 1. Plot of surface density S (in code units) against radius for a spreading Gaussian disc with: (a) no advective term and; (b) a 
strong advective term which results from a magnetic field. The evolution of the depleted region or 'hole' is clearly visible in the latter 
case. 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 8.5 9 9.5 10 

R/Rco R/Rco 



Figure 2. Plot of mass transfer rate, M (in code units), against radius for a spreading Gaussian disc with: (a) no advective term and; 
(b) a strong advective term which results from a magnetic field. Where £ < 0, M is not plotted. These are the same simulations, plotted 
at the same times as Figure HI 



graphs to show the position in the disc where surface density 
first exceeds S^f t h and outbursts would begin. 

Figure a) represents a simulation of the disc with no 
magnetic field. The inner radius of the simulation lies at the 
surface of the star R* ~ 1 R©. Plot (a) is on the same scale 
as plots (b) and (c) for ease of comparison, and plot (a)i is a 
zoomed in view of graph (a). It can be seen in plot (a) that 
a fairly flat surface density profile is formed with an upturn 
in surface density developing towards the surface of the star. 
The surface density rises with time and the critical density 
is first reached at the surface of the star, as expected. 

Also in Figure the effects of weaker and stronger pro- 
pellers can be seen in plots (b) and (c) respectively. In both 
cases the early stages of disc evolution are not dissimilar to 
those in case (a) for most of the disc. The surface density 
increases with time and as it does so the effect of the mag- 
netic field becomes apparent. In case (b) the inflow of mass 



towards the star is suppressed inside ~ 272 co and the critical 
density is first exceeded at ~ 4ii co . This occurs much later 
than in plot (a), as expected. 

In plot[3{c) a stronger propeller is illustrated, where the 
advection term is greater. The plot is qualitatively very sim- 
ilar to plot (b) but the completely depleted region extends 
to 3i?co and the outburst would begin at ~ 5R co ~ 50 Rq. 
Although a full outburst cycle is not simulated and therefore 
it is not possible to draw detailed conclusions about i rcc , it 
can be seen that outbursts would begin later with increasing 
P in these simulations. Values of R\ of the expected order 
of magnitude can be seen in Figure 03 If the out bursting re- 
gion has a width of order a tenth of the radius at which the 
outburst begins, as assumed in Section [5] then the approxi- 
mation that E ~ E cr it in this region is clearly a good one in 
all three cases. 

The mass transfer M in the disc for all these simulations 
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Figure 3. Plots of surface density S (in code units) against radius where advective terms of varying strengths, representing magnetic 
fields, have been applied. In all cases the initial conditions were an empty disc with a constant inflow of M ~ 10 -7 M0yr— 1 from the 
right. Plot (a) shows the behaviour of a disc with a zero magnetic field while plots (b) and (c) represent a weak and a stronger field 
respectively. In both magnetic cases an inner disc truncation can clearly be seen, although this occurs at a greater radius in case (c). 
Plot (a)i is a 'zoomed in' version of plot (a) showing where the density first exceeds £j^? t h , while plots (a), (b) and (c) are shown on 
identical axes for ease of comparison. In each case the latest plot is at the point where E first reaches the critical density. 



is plotted in Figure0]at the same times as in Figure[3] In all 

cases the disc has yet to reach steady state at the point where 
the outburst begins. If no outburst occurred, the disc would 
eventually approach steady state. Because the spreading of 
the disc slows as it nears steady state, the longest recurrence 
times will be observed when the critical density is reached at, 
or shortly before, steady state. When (3 is increased so that 
steady state is reached with the surface density maximum 
at larger radii then no outbursts will occur and the disc will 
remain quiescent indefinitely. Such a steady state without 
outbursts could exist for example, in our results, when (5 > 
1 x 10~ 13 . This is equivalent to a stellar magnetic field of 
B > 1 kG in a typical CTT disc with a c ~ 0.01. 

It is possible to observe trends in both Ri and t lac as 
a function of j3. Several simulations were performed with 
different values of the parameter j3. As expected, the inner 
radius increases with j3. This increase can be fitted very well 
with a power law of index ~ 1/3. The results are shown in 



Figure The values of Ri obtained here are consistent with 
those which are required in Section |5] 

In Figure 6 we plot the relation between j3 and time un- 
til first outburst, starting with an empty disc. In plot (a) the 
time before outburst increases rapidly beyond a certain f3. 
This is a result of the disc model approaching steady state. 
This means that it is possible to obtain an arbitrarily long 
recurrence time provided the star has a sufficiently strong 
magnetic moment. The value of f3 at which this occurs is a 
function of 8 but the form remains the same. Plot (b) shows 
an enlarged section of plot (a). We additionally plot the cor- 
responding inner disc radius. The relation between Ri and 
the time before first outburst can be fitted by an expression 
of the form t oc Rf as would be expected from Equation 
(5). Further, when the time before outburst is long, the in- 
ner radius of the disc is close to Ri ~ 40 Kq as predicted 
in Section 2. However these results do not provide quanti- 
tative estimates of t lac because the initial conditions are so 
different from those just after an outburst. 
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Figure 4. Plots of M (in code units) against radius for the same three advective terms strengths as were used in Figure 13 The input M 
is plotted for comparison. Once again, where S < 0, M is not plotted. If no outburst occurred then the disc would eventually approach 
steady state. This can be seen in the plots at 500 yr, where the disc has been allowed to evolve without outburst. 




4e-14 6e-14 
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Figure 5. Plot of inner disc radius i*?; against magnetic parameter /3 for a viscosity parameter of S = 2 X 10~ 8 and /3 in code units. The 
curve can be fitted well by a power law of the form Ri oc /3 1 / 3 . 
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Figure 6. Plots of time before first outburst against magnetic parameter (3. In plot (a) the time before the outburst increases very rapidly 
beyond a certain j3. There, steady state is nearly reached before outburst. Plot (b) is a 'zoomed in' view of plot (a), also indicating the 
corresponding inner disc radii. 
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Figure 7. Plot (a) shows temperature in Kelvin as a function of radius in steady state disc models with and without magnetic truncation. 
The hottest part of the disc is missing in the truncated case. Plot (b) shows the spectral energy distribution in arbitrary units against 
wavelength in microns. In the truncated case the flux is reduced for wavelengths below 10/xm. 



4.3 Simulated spectra 

In order to compare the results of numerical models to ob- 
servations it is useful to produce simulated spectra. This can 
be done in a simple way by dividing the disc in to a series of 
concentric annuli. The viscou s dissipation in eac h annulus 
can then be approximated by l|Frank et al. 1 12002*) 



E? 



(29) 



where v is the Shakura-Sunyaev viscosity as given by Equa- 
tion l|24[l and E is taken from simulation results. Stellar irra- 
diation may effect the dis c temperature and he nce the SED. 
Following the approach of lAdams fc Shul il986T) the total lu- 
minosity of the disc can be broken down into components 
due to viscous dissipation and to disc ir radiation. The vis - 
cous component can be approximated by jFrank et al. .120021) 



S v = 



GM+M 
2R+ 



(30) 



In the limit of a large disc, the component of disc luminosity 
due t o irradiation can be approximated by l|Adams fc Stml 
1986) 



(31) 



where T* denotes the temperature of the star. The two 
components can be compared using typical CTT parame- 
ters of T* ~ 4000 K, M* ~ 1 M , M ~ 10" 7 Myr -1 and 
R* ~ 1 R,0. In this case, the component of disc luminosity 
due to irradiation is of the order of a few percent of that due 
to viscosity. It is therefore resonable to neglect irradiation in 
the simple model applied here. Considering viscous heating 
only and assuming that the disc is in thermal equilibrium, 
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the temperature of each annulus can be calculated from the 
Stefan-Boltzmann law according to 



D (R) = eaT 4 



(32) 



where e is emissivity (assumed to be unity), a is the Stefan- 
Boltzmann constant and T represents the effective temper- 
ature. Temperature profiles, calculated using this simple 
model, can be seen in Figure Efa) for both the truncated 
and non-truncated disc models. A spectrum emitted from 
each annulus is then generated, assuming that the disc be- 
haves as a black body. The spectral radiancy is then given 
by the Planck radiation law so that 



R v — 



2nhv A 



(33) 



where h represents Planck's constant, v is the frequency of 
the radiation, k is Boltzmann's constant and c is the speed 
of light. All the spectra of the annuli are summed to give 
the spectral energy distribution (SED) of the entire disc. A 
result is presented Figure [7Jb). In the case of the truncated 
disc the flux at wavelengths shorter than lfim is greatly re- 
duced. This is to be expected, as the hottest part of the disc 
has been significantly depleted. This may help to explain 
the infrared excess observed in some YSOs. The shape of 
the SED at longer wavelengths is not fully resolved as the 
contributions of the outer disc, beyond 200 R.0, have not 
been assessed. 

Observed spectral energy distributions of CTTs tend 
to include an infrared excess which is attributable, at least 
in part, to the presence of a circumstellar accretion disc 
iBertout . Basri fc Bouviei]ll988t lAdams. Emerson fc Fuller! 
Il99d) . Iltoh et al.l fcOolT successfully reproduce observed 
SEDs for CTTs by summing the SED expected from the 
young stars themselves with disc SEDs. The position in 
wavelen gth of the inner e dges, and peaks, of the disc SEDs 
M 



used bv lltoh et al.l (|2003j) lie between those of our truncated 
and non-truncated models which have been al lowed to reach 
the st eady state without reaching outburst. lUltchin et al.l 
( 2001) argue that a truncated accretion disc produces an 
infrared excess more closely in agreement with observations 
than an unmodified disc model. It should be noted that it is 
very difficult to separate the disc component from the domi- 
nant stellar component of the observed SEDs at wavelengths 
shorter than 10/im. Detailed comparison with observation is 
problematic as the method used to generate disc SEDs from 
density profiles in this paper is only very approximate, and 
takes no account of material which is able to reach radii less 
than J? mag . Further comparison of observational results to 
the predictions of our model would therefore require a more 
detailed analysis than is possible under the restrictions of a 
one dimensional treatment. 



5 DISCUSSION 

The FU Ori stars undergo what appear to be thermal- 
viscous outbursts. It is very difficult however to reconcile 
observational estimates of t lec and t r i sc with the very short 
estimates obtained from the standard thermal-viscous disc 
instability model with a c > 10~ 3 It is possible to do so 
by using a very low value of a c . We propose an alternative 



scheme in which outbursts occurring in an accretion disc 
truncated at a radius Ri ~ 40 R.0 are consistent with the 
observational estimates of the outburst, recurrence, and rise 
times. The central star, acting as a magnetic propeller, can 
produce the necessary truncation of the disc. The results 
of our one dimensional model confirm the feasibility of this 
idea. A depleted region in the centre of the disc can be made 
large enough to produce t rcc similar to those estimated from 
observations. In particular, an arbitrarily long t rcc can be 
obtained by selecting appropriate values for viscosity and 
magnetic field strength. 

The simulation of a full outburst in a refinement of 
the one dimensional code described above, or using a three 
dimensional technique such as smoothed particle hydrody- 
namics, would provide more accurate estimates of the out- 
burst timescales, which could be compared to analytical 
results and to observations. Three dimensional simulations 
would be particularly valuable as they would produce a self 
consistent picture of the inner disc. It would also be desirable 
to treat viscosity in a more rigorous manner. A more exten- 
sive study of stellar spin evolution would allow quantitative 
analysis of any possible feedback mechanisms. 

The prolongation of recurrence times in accreting ob- 
jects as a result of the magnetic propeller mechanism need 
not be confined to young stellar objects. The same mecha- 
nism could potentially exist in some accreting binaries. In 
YSOs the full significance of the magnetic propeller may not 
yet be known. It is possible for example that a depleted inner 
disc may have an effect on the planet formation process. 
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